Threshold hyperon production in proton— proton collisions at COSY— 11 
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E + hyperon production was measured at the COSY-11 spectrometer via the pp — ^ nK + Y, + reac- 
tion at excess energies of Q = 13 MeV and Q = 60 MeV. These measurements continue systematic 
hyperon production studies via the pp — > pK + A / E° reactions where a strong decrease of the cross 
section ratio near threshold was observed. In order to verify models developed for the description of 
the A and E° production we have performed the measurement on the E + hyperon and found unex- 
pectedly that the total cross section is by more than one order of magnitude larger than predicted 
by all models investigated. 

After the reconstruction of the kaon and neutron four momenta, the E + is identified via the 
missing mass technique. Details of the method and the measurement will be given and discussed in 
view of theoretical models. 

PACS numbers: 13.60.Hb, 13.75.-n, 25.40.Ve, 28.20.-n 
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I. INTRODUCTION 



The study of hyperon production in hadron induced 
multi particle exit channels like pp — > NKY involves 
several aspects. The nucleon-hyperon interaction can be 
extracted by analyzing the NY subsystem in the appro- 
priate kinematical region. Closely related to that is the 
issue of the reaction mechanisms of the hyperon produc- 
tion which have to be clarified for an unambiguous inter- 
pretation of the data. If the hyperon production is due 
to the excitation and a subsequent decay of intermediate 
nucleon resonances than information about the structure 
of the relevant resonances can be extracted. 

The pp — > pK + A excitation function near threshold 
shows a clear deviation from a pure s-wave phase space 
distribution, and the proton-hyperon final state interac- 
tion (FSI) has to be included to describe the data [HQ-El 
4]. In the pp — > pK + TP channel the pY FSI seems to be 
negligible and the pure s-wave phase space calculations 
follow the data points reasonably well. The cross section 
ratio a(pp — > pK + A )/a(pp — > pK + TP ) below excess en- 
ergies of Q~ 20 MeV is about 28 0,0] in contrast to the 
value of about 2.5 determined for excess energies higher 
than Q = 300 MeV [j| (see figure^) . This value is in good 
agreement with the A/E° isospin relation. The question 
arises if this drastic increase of the cross section ratio 
near threshold is a mere effect of the NA FSI or wether 
it is partly due to the reaction mechanisms in the NY 
channels. To explain the increase of the A/E° cross sec- 
tion ratio in the near threshold region, different scenarios 
were proposed. 
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FIG. 1: The cross section ratio for E° and A production in 
the threshold region. The data Q, 0, 0, @ are compared to 
different model predictions. The curves are theoretical pre- 
dictions described in the text. 



Calculations of the strangeness production based on 
incoherent it and K exchange have been performed in 
Ref. @. The irN -> YK and KN -> YK scatter- 
ing amplitudes for pion and kaon exchange, respectively, 
were taken from the existing data in the higher en- 
ergy region [s]. This incoherent tt/K exchange model 
describes the pp — > pK + A cross section over the whole 
energy range, however in the near threshold region the 
pp — > pK + Y,° channel is overestimated, and thus the pre- 
dicted A/S° ratio is too low for Q ^ 20 MeV (see curve 
(1) in figure nj. 

A better description of the strong rise of the ratio to- 
wards lower Q- values is achieved by the resonance model 
(curve (2) in figure [TJ) |E SB l3- In this model the 
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nonresonant direct contributions like ir or K exchange 
were not included, but the n, 77 and p meson exchange 
with the excitation of the intermediate baryonic reso- 
nances JV(1650), #(1710), #(1720) and A(1920) are 
taken into account. In this resonance model the near 
threshold region of the A/E° cross section ratio seems 
to be better reproduced than the higher energy values 
(i.e. Q 10 MeV). It should be stressed that in these 
calculations the parameters were fixed on the basis of 
higher energy data, before the near threshold A and E° 
data were available. 

Other calculations by Shyam ^(j (based on the ef- 
fective Lagrangian model) for the strangeness produc- 
tion include meson exchange (tt, p,o~ and u>) together 
with the excitation of resonances. The coupling con- 
stant was determined from data of the n + p — > E + #+, 
ir~p — > E°#° and ir~p — > E~ K + reaction channels. The 
coherent sum of resonant states and meson exchange pro- 
cesses describe the experimental data for the pp — > pK + A 
and pp — > pK + YP channels very well. The effective La- 
grangian model is depicted by the curve (3) in figure ^ 

The Jiilich theory group has performed calculations in- 
cluding 7r and K exchange 0,^2- -"- n their approach the 
interaction between the hyperons (A, E) and the nucleon 
is described by a microscopic (AN-EN) coupled chan- 
nel model |13| with a coherent superposition of the pro- 
duction amplitudes. The A production is dominated by 
the K exchange and therefore the contribution due to 
an interference between n and K exchange is negligible 
in this hyperon channel. On the other hand the 7r and 
K exchanges give a comparable contribution to the cross 
section in the case of E° production. An interference be- 
tween 7r and K exchange amplitudes act very differently 
on the two channels. Within the Jiilich meson exchange 
model the large cross section ratio can be described by a 
destructive interference of the 7r and K exchange ampli- 
tudes only. For excess energies above 20 MeV the model 
is not valid any more but qualitatively the cross section 
ratio given by the model stays at a nearly constant level. 

Although the various descriptions of the cross section 
ratio differ even in the dominant reaction mechanism, 
all reproduce more or less the trend of an increase of 
the A/E° cross section ratio in the threshold region (see 
figure The present data are not sufficient to un- 
ambiguously identify the dominant reaction mechanism. 
To clarify this point further data are needed. Espe- 
cially the other isospin channels should allow informa- 
tion about the production mechanisms to be extracted. 
Recently, after the installation of a neutron detector in 
addition to the A and E° production the reaction chan- 
nel pp — > n# + E + became accessible at the COSY-11 de- 
tection system. The measurement of the E + hyperon 
production via this reaction was performed at two beam 
momenta, Pbeam = 2.6 GeV/c and Pbeam = 2.74 GeV/c, 
corresponding to excess energies of 13 MeV and 60 MeV, 
respectively. 



II. EXPERIMENT 

COSY-11 is an internal magnetic spectrometer at the 
COoler SYnchroton and storage ring COSY [14| in Jiilich. 
The interaction between a proton in the beam and a pro- 
ton from the H2 cluster target may lead to the pro- 
duction of the E + hyperon in the pp — * n# + E + reac- 
tion. The charged reaction products are separated from 
the circulating beam in the magnetic field of one of the 
regular COSY dipoles ^(|. The E + hyperon is identi- 
fied via the missing mass technique by detecting the K + 
and the neutron. The momentum vector (p) of the K + 
meson can be determined by tracking the K + trajectory 
reconstructed in the drift chambers (DC1 and DC2 in 
figure |2l through the known magnetic field back to the 
target point. The momentum vector together with the 
velocity (/?) measurement in the two scintillators S8 and 
SI are used to identify the kaon via its invariant mass 

(Mf nv =^(/r 2 -i)). 

Assuming that a hit in the neutron detector is from 
a neutron, the four momentum vector of the neutron is 
calculated by the measured velocity, the direction of the 
neutron (given by the first module hit) and the known 
mass. The background from charged particles incident 
on the neutron detector is discriminated by veto scintil- 
lators. 
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FIG. 2: COSY-11 detection set-up [TtJ with the superim- 
posed tracks of kaons and neutrons from the pp — > nK + E + . 
The scales of the axes are given in [cm]. 



In figure the experimental distributions of the 
squared missing mass (m~.) of the pp — ► nK + X system 
for the two beam momenta are shown. For the higher mo- 
mentum, an enhancement around the squared E + mass 
is clearly seen on a large background (figure |2t>), but for 
the lower beam momentum (figure [3^,) a E + peak is not 
directly visible. 
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FIG. 3: The experimental squared missing mass spectra of 
the pp — » nK + X system for two investigated beam momenta. 
The arrows point to the nominal squared mass of the E + hy- 
peron. The statistical errors of the counting rates are shown. 



pp — > pK + A"f and pp — > pK + A"/~/ are the dominant back- 
ground channels in the S + region. The Monte Carlo code 
includes the realistic geometry and physics processes like 
energy loss and straggling which occasionally cause the 
misidentification of the particle type. 



All background channels result in a rather smooth dis- 
tribution of the missing mass spectrum as can be inferred 
from calculations from Monte-Carlo studies and by com- 
paring the two experimental distributions (see figure 



In order to determine the number of S + events in 
the higher energy data set, the background has been 
parametrized by a polynomial function which is added to 
the expected missing mass distribution of the nK + sys- 
tem for the pp — > nK + T, + reaction obtained from sim- 
ulation studies. In figure 0^ the experimental missing 
mass spectrum of the pp —> nK + X system is compared 
with the fitted polynomial function. The expected dis- 
tribution from Monte-Carlo studies with X = S + is de- 
picted in the figure as well. Figure^ shows the result of 
the subtraction of the polynomial from the experimental 
missing mass distribution together with the Monte-Carlo 
distribution. 



For the lower energy data set a S + peak is not obvi- 
ously visible via the missing mass distribution. There- 
fore, a simple polynomial background parametrization 
cannot be used. To determine the number of S + events it 
was assumed that the background shape for this data set 
is the same as that at the higher energy. This assumption 
is justified since there are no new open channels for the 
higher energy. 



(a) (b) 
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FIG. 4: Background determination for the pp — » nK + T, + re- 
action at Pbeam = 2.74 GeV/c. (a) The experimental squared 
missing mass spectrum of nK + system with a polynomial 
background fit and the simulated pp — > nK + T, + spectrum, 
(b) Result after subtracting the background from the experi- 
mental distribution in comparison to the simulated spectrum. 
The arrows show the nominal squared mass of the S + hy- 
peron. Only statistical errors are shown. 

In order to understand the background distribution, 22 
reaction channels (mostly multi-pion reactions but also 
pp — » pK + A (E°, A7)) were simulated and their contribu- 
tions to the missing mass distribution were determined. 
These studies showed that the reactions pp — ► pK + A , 



At the COSY-11 experiment the shape of the missing 
mass distribution is mainly determined by the acceptance 
of the detection system and is dependent on the excess 
energy of an individual event. From the analysis of r\ 
and 77' production studies at COSY-11 it was verified 
that the background shape resulting here mainly from 
multi pion production is in very good agreement with 
the expectations from Monte Carlo studies taking into 
account the detector characteristics and is comparable at 
different beam momenta. In addition, Monte Carlo data 
of the reaction channels which contribute dominantly to 
the background in the S + production were compared in 
view of the background shape by adjusting the kincmat- 
ical limits. Within error bars their shapes were identical. 
Therefore it is justified to assume that the background 
shape is the same for both beam momenta. For a de- 
tailed discussion on the background shape at COSY-11 
we refer to 17]. The background shape from the exper- 
imental missing mass distribution for the higher energy 
data set was shifted such that the kinematical limits for 
both spectra were the same. Afterwards, it was added to 
the missing mass distribution of the pp — ► nK + H + reac- 
tion which was determined from Monte-Carlo studies at 
the lower energy (see figure |SJ). The relative magnitudes 
were determined by the \ 2 minimisation of the functional 
forms to the data sample at lower energy. 
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FIG. 5: Background investigation for the pp — > nK + 'E + re- 
action at Pbeam = 2.6 GeV/c. (a) The experimental squared 
missing mass spectrum of the pp — » n7\ + X system is com- 
pared to the assumed background taken from the data at Q 
= 60 MeV and the simulated pp nK + T, + distribution, (b) 
Result of the subtraction of the background from the experi- 
mental distribution in comparison to the simulated spectrum. 
The arrows show the nominal squared mass of the S + hy- 
peron. Only statistical errors are shown. 



III. RESULTS 

A. Total cross section 

For the lower energy data set, no clear enhancement 
near the kaon mass in the invariant mass distribution is 
visible. Therefore the assumption was made, that the 
kaon peak in the experimental distribution has the same 
position and width as in the simulated distribution. 

As a cross check, event samples in different regions of 
experimental invariant mass but still within the expected 
range for kaons were taken and the missing mass distribu- 
tions were determined. The assumed background shape 
was subtracted from the experimental distribution and 
the enhancement around the mass of the S + was inter- 
preted as a signal from the S + hyperon (see figure EJ). 
This procedure was repeated for five different invariant 
mass regions and the number of events under the remain- 
ing peak was determined. The results are shown in fig- 
ure |SJ The vertical error bars correspond to the statis- 
tical error of the number of events. The horizontal er- 
ror bars indicate the widths of the appropriate invariant 
mass intervals. In figure [B] curve (1) shows the expected 
invariant mass distribution of kaons from Monte-Carlo 
calculations. Curve (2) is a Gauss fit to the experimen- 
tally determined number of £ + events in the different 
invariant mass regions. This can be interpreted as an 
experimental distribution of kaons. The assumed experi- 
mental kaon invariant mass distribution is consistent with 
the distribution expected from Monte-Carlo studies. 
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FIG. 6: The number of E events obtained for different 
intervals of the K + invariant mass. Curve (1) - expected 
from the Monte-Carlo studies distribution of the K + events. 
Curve (2) - Gauss fit to the points on the figure given by the 
number of S + events. See text for details. 

For the higher beam momentum a kaon peak is clearly 
visible in the invariant mass distribution and can be fit 
by a Gauss function. In a systematical study of the result 
at the higher beam momentum three different selections 
on the experimental invariant mass distribution were ap- 
plied. Events from the following regions: /i exp i0.25a exp , 
[i exp ± 0.5a exp and fi exp ± 1.0a exp centered around the 
nominal reconstructed mass (fi e xp) were taken and the 
corresponding missing mass distributions were generated. 
a exp is the experimental resolution of this observable. 
Next the number of S + events for each of these distribu- 
tions was determined. The results are listed in the middle 
column of tableland in the last column the number of 
the £+ events corresponding to the full Gaussian distri- 
bution. 



TABLE I: The number of E + events as a function of the 
width of the invariant mass range. Only statistical errors are 
given. 



Invariant 
mass range 

f-l>exp 


Number 
of events 
identified as E + 


Number of events identified 
as E + and scaled to the 
full Gaussian distribution 


± 0.25 a exp 


187 ± 51 


944 ± 257 


i 0.5 (J exp 


367 ± 83 


960 ± 217 


± 1 .0 (Jexp 


661 ± 134 


969 ± 197 



In order to calculate the cross section for the 
pp — > nK + T, + reaction the number of E + hyperon events 
and the detection efficiency of the COSY-1I apparatus 
for the two excess energies were determined. The lumi- 
nosity was determined by a simultaneous measurement 
of proton-proton elastic scattering. 

In table ITT1 the total cross sections for both beam mo- 
menta are given. The systematical errors are due to: i) 
error of the detection efficiency determination which is 
8.5 % for the lower and 3.5 % for the higher energy data 
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set (including the inaccuracy of the effective detector po- 
sition and of the beam momentum determination), ii) 
uncertainty in the form of the background, and Hi) er- 
ror of the luminosity calculation which is 3 % for both 
data sets and includes the uncertainty due to the nor- 
malization procedure and the error of the solid angle 
determination. For the data at 2.74 GeV/c the uncer- 
tainty in the background form was estimated by com- 
paring the polynomial parametrization to a background 
form resulting from the sum of known background re- 
action channels generated in Monte Carlo studies. The 
difference is about 18%. For the data at 2.6 GeV/c the 
m\ region used to normalize the background was varied 
resulting in an error of about 20%. The values in table ITT1 
include also a change in the detection efficiency resulting 
from the inclusion of higher partial waves. Close to the 
reaction threshold partial wave contributions higher than 
L = are not expected. However, if the excess energy 
Q for a given channel goes beyond a few MeV it cannot 
be a priori neglected. The contribution of partial waves 
with L > is not known for the reaction investigated 
here. Thus, this effect has been studied by assuming a 
similar distribution of partial waves as measured in the 
pp — > pK + A channel studied by the TOF collaboration 
at COSY jH [H HJ. The inclusion of higher partial 
waves with a relative strength quoted in 0, ^| lower 
the detection efficiency by 30 % for the lower and by 7.7 
% for the higher energy data set. The sum of the total 
systematical errors equals 60 % for the lower and 34 % 
for the higher energy data set. 

TABLE II: Total cross section for the reaction pp — > nK + Y, + . 
Statistical and systematical errors are presented, respectively. 



Beam momentum 

Pbeam 

[GeV/c] 


Excess energy 

Q 

[MeV] 


Total cross section 
a 

[pb] 


2.60 


13 


4.56 ±0.94 ±2.7 


2.74 


60 


44.8 ± 10.7 ± 15.2 



B. Comparison with model predictions 

Among the models described in the introduction only 
two give predictions for the pp — > nK + T, + reaction, 
namely the Jiilich meson exchange model [Til . Il2j | and 
the resonance model 0, . Calculations of the S + pro- 
duction within the Jiilich meson exchange model predict 
a total cross section of a — 0.23 fib at Q = 13 MeV for 
the destructive interference (which was necessary to de- 
scribe the high A/E° cross section ratio at threshold). 
This is about a factor of 20 below the experimental value 
of 4.56 pb given in table ITT1 A constructive interference 
would result in a cross section even a factor of 53 too low. 

For the resonance model the predictions for the 
pp — > nK + T, + channel for the near threshold region de- 
viate even more from the data. In figure {7\ the model 
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FIG. 7: Comparison of the experimental total cross sec- 
tion with the resonance model [(j, 0J predictions for various 
pp — » NK + Y reactions. Full triangles in (a) are data obtained 
in this work. Data in the near threshold region (presented as 
full symbols in (b) and (c)) are taken from Refs. Qlftu and 
data from the high excess energy region (open symbols) from 
Refs. 0, 21]. In (c) the data point indicated by the arrow was 
determined from our data as a cross check of the luminosity 
calculation. 



predictions and the available data for the pp — > nK + T, + 
(a), pp — ► pK + YP (b) and pp — > pK + A (c) channels are 
shown. The filled data points presented by triangles, dots 
and squares in the near threshold region were measured 
by the COSY-11 collaboration dUQ. The data point 
in the pp — > pK + A channel indicated by the arrow was 
determined in parallel to this work by selecting the K + p 
exit channel which was included in the events sample at 
2.74 GeV/c. The high energy data for the given reactions 
were taken from |5J and [2l| . The model calculations for 
each channel are given by the solid lines 0, • 

The data point for the pp — » nK + Y, + channel at 
Q = 13 MeV is underestimated in the total cross sec- 
tion calculated using the resonance model [6J, |9j by about 
a factor of 500 and for Q — 60 MeV by about a factor 
of 50. For the pp — ► pK + TP channel, this model calcu- 
lation describes the existing data set and in the case of 
the pp — > pK + A channel the underestimation of the cross 
section in the near threshold region is about a factor of 
16 being 30 times smaller than for the E + production. At 
high excess energies, the S + data points are by factor 3 - 
4 below the model calculations. Previous COSY-11 hy- 
peron production studies conclude that final state inter- 
actions (FSI) plays an important role in the near thresh- 
old A production |lll2j.l3 Ll23l |. In the resonance model the 
FSI is not included |6|, |9( and therefore the deviation of 
the model calculations from the data points in the near 
threshold region is expected if a strong FSI is present. 
This effect is clearly seen for the pp — > pK + A and barely 
observed for the pp — > pK + YP reaction channel. 
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FIG. 8: The pp -> nA"+S+ , pp -> pK + A and pp -> pif+E° 
cross sections as a function of the excess energy Q. Experi- 
mental data are from Refs. Q, 0, S El and from this work. 
The error bars for the pp — > nK + T, + data represent a sum 
of statistical and systematical uncertainties given in table ITT1 
The lines show the calculations corresponding to 3-body phase 
space with (solid line) and without (dashed line) final state 
interaction. 

Investigations of hyperon production in COSY-11 
|H S Q indicate that the pp — > pK + A data cannot be 
described by a pure 3-body phase space (PS) dependent 
cross section expressed as [25j : 

a = K-Q 2 , (1) 

where K is a normalization factor and Q the excess en- 
ergy. Therefore a modification is needed which takes into 
account the proton - hyperon FSI. In order to describe 
the near threshold region, the parametrisation of the ex- 
citation function including the FSI proposed by Faldt- 
Wilkin U was used. It is expressed by: 



(i + v /iTq7^ 2 ' 

where C and e are parameters related to the FSI strength. 

In figure |H1 the cross sections for different production 
channels for the A, S° and S + hyperons are compared to 
predictions based on the 3-body phase space (PS, dotted 
line) and the 3-body phase space calculations modified 
by the pY FSI (PS + FSI, solid line), following equation^ 
with e and C as free parameters. These parameters are 
related to the scattering length a and the effective range 
r of the pY potential Q . 

For the pp — > nK + Y, + data the resulting e and C 
parameters have similar values as for the pp — > pK + A 
channel. It seems that for the S + production via the 
pp — > nK + H + reaction a rather strong n — S + FSI is 
present, however, within the error bars also the curve ob- 
tained without FSI describes the two data points. There- 
fore an unambiguous conclusion about the n — E + FSI 



requires more data to disentangle the reaction mecha- 
nisms and especially the role of nuclear resonances. 



IV. CONCLUSIONS AND PERSPECTIVES 

The total cross section of the pp — > nK + Y, + reaction 
was determined by means of the COSY-11 detection sys- 
tem for excess energies of Q = 13 and 60 MeV. The values 
established are more than an order of magnitude larger 
than the expectations of currently available model pre- 
dictions. 

The unexpected large total S + production cross sec- 
tion is consistent with an observation by Tan |26j| who 
concluded that when assuming charge symmetry in S + n 
and S° p scattering, the contribution from the E° dia- 
gram is less than one seventh of the contribution from 
the S + channel. Further, recently it was suggested 
that for the <fi production [27| a strong enhancement of 
the reaction amplitude towards threshold might be due 
to the presence of a crypto exotic baryon with hidden 
strangeness. Although this observation is not directly ap- 
plicable to other isospin channels, (in the T,°p system no 
corresponding structure was observed) it might be a hint 
of an exotic mechanism. In any case, present theoretical 
predictions of the cross sections strongly underestimate 
the experimental data. The comparison of the excitation 
function expected based upon phase space distribution 
including N — Y FSI to the data results in FSI parame- 
ters comparable to the p — A system. This may indicate 
a strong n — S + interaction but due to the large system- 
atic uncertainties the data are also consistent with a pure 
phase space distribution without n — S + FSI. 

Further studies of the E + production are necessary to 
clarify the picture. On the experimental side additional 
data points should be added for which an improved event 
selectivity is favorable to reduce the large uncertainties 
introduced by the background subtraction. A 4-7T de- 
tection system for neutral and charged particles which 
will be soon available with WASA at COSY could be 
used |28f. On the theoretical side an improved model 
has to be developed which consistently reproduces the 
hyperon cross section data near threshold. 
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VI. NOTE ADDED IN PROOF 

During the evaluation process of the article we have 
been made aware of the predictions of the excitation func- 
tion of the total cross section for the pp — > nK + Y, + re- 



action which is closer to the data in comparison with 
the models discussed, yet still underpredict the deter- 
mined total cross sections by more than an order of mag- 
nitude 0. 
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